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Abstract 



Using a description of defects in solids in terms of three-dimensional gravity, 



H ' we study the propagation of electrons in the background of disclinations and 

screw dislocations. We study the situations where there are bound states that 
are effectively localized on the defect and hence can be described in terms of 
an effective 1 + 1 dimensional held theory for the low energy excitations. In 
the case of screw dislocations, we hnd that these excitations are chiral and 
can be described by an effective field theory of chiral fermions. Fermions of 
both chirality occur even for a given direction of the magnetic field. The "net" 
chirality of the system however is not always the same for a given direction of 
the magnetic field, but changes from one sign of the chirality through zero to 
the other sign as the Fermi momentum or the magnitude of the magnetic flux 

1 



is varied. On coupling to an external electromagnetic field, the latter becomes 
anomalous, and predicts novel conduction properties for these materials. 

I. INTRODUCTION 

The effect of defects on properties of solids has been a very active field of research . The 
fact that the theory of defects in solids can be reformulated as a version of three-dimensional 
gravity has been discussed in ||^. This formulation corresponds to looking at the contin- 
uum limit of a crystalline solid in which the static defect configuration is characterized by 
a non-trivial metric corresponding to a static spacetime. There is a one-to-one correspon- 
dence between the classification of defects in terms of Burgers vectors and different metric 
configurations. This formulation has the merit of highlighting the geometrical properties of 
defect configurations. 

The properties of electrons interacting with these static defects can be studied by looking 
at their propagation in the background of these defects 0]. In this paper, we will focus on 
situations where one can have bound states of electrons localized on line defects like the 
disclination and the screw dislocation. In this context, it is useful to recall that there is a 
very close similarity of this system to that of particles in the background of cosmic strings and 
domain walls . It is known that there exist states localized on strings and domain walls, 
leading to interesting phenomena like gauge and gravitational anomalies. The properties of 
states localized on domain walls in the context of condensed matter systems have also been 
studied in [^,0]. 

As mentioned above, we here study situations where there are bound states localized 
on line defects. In cases where one has chiral bound states, the effective description of low 
lying excitations at the defects is given in terms of a (1 + 1) dimensional massless chiral 
Dirac fermion theory. Using well-known results about the gauge theory of chiral fermions 
in (1 + 1) dimensions, one can show that there is a U{1) anomaly [§, leading to charge 
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nonconservation in the (1 + 1) dimensional effective action. One can also show that the 
charge nonconservation is exactly cancelled by the neglected bulk modes These modes 
contribute to the bulk current. The latter leads to a net inffow of charge to the defect that 
precisely accounts for the charge violation given by the U{1) anomaly. This effect of the 
bulk modes can be summarized in terms of effective Chern-Simons actions on planes having 
the defect as boundary. 

In order to give a simple but explicit example of the above scenario, we first study the 
Schrodinger equation of the electron in an external electromagnetic field in the absence 
of any defect. The coupling of the latter to the magnetic moment of the electron is also 
taken into account. We assume that the gyromagnetic ratio of the electron is two. In 
this case, it is known that the Dirac Hamiltonian in two dimensions can be squared to 
obtain the Schrodinger equation with magnetic moment interaction. The corresponding two 
dimensional (2d) Schrodinger Hamiltonian is the part of the 3d Schrodinger Hamiltonian 
transverse to the defect. We first analyze the 2d Dirac Hamiltonian and its domains of self- 
adjointness. The latter are found to be compatible with the domains for the 2d Schrodinger 
Hamiltonian when the derivatives of the Dirac wave function are in the same domain as the 
wave function itself. 

In the absence of defects, the above 2d Dirac equation and the corresponding 2d 
Schrodinger equation are known to have bound states of zero energy. The number of bound 
states is given by an index theorem and is related to the integer part of the flux through the 



plane. This is the Aharonov-Casher theorem [jTO|. These bound states are localized in the 
region of the magnetic field. 

We extend the Aharonov-Casher analysis to situations where there are special types of 
defects, namely the edge disclination and the screw dislocation. We find that similar results 
exist in this situation, that is we can extend the counting of bound states to the situation 
where there is a conical deficit angle and a screw defect. 

Next we obtain a (1 + 1) dimensional field theory for the bound states as they move along 
the defect by linearizing the excitations about the Fermi surface. It is a (1 + 1) -dimensional 



chiral fermion theory for the low energy excitations. When this effective field theory is 
coupled to an external electromagnetic field, there is then an anomaly |p . We then proceed 
to look at the response of the system to a constant electric field applied in the direction of 
the defect and find that there is a net inflow of charge from the planes transverse to the 
defect, in other words a "transverse Hall effect", because of the anomaly. 

The paper is arranged as follows. In section 2, we discuss the problem of the motion 
of an electron on a plane (with or without a point removed ) in an external magnetic field 
and find the Dirac Hamiltonian. In section 3, we do the partial-wave analysis for this 
Hamiltonian when the plane has a hole and find its domains of self-adjointness. We also find 
the zero-energy solutions for the Dirac Hamiltonian. We introduce the line defects, namely 
disclinations and screw defects, in section 4. We also find the zero energy solutions in these 
background geometries, which are very similar to their flat space counterparts. These states 
turn out to be individually chiral in terms of their motion around the defect, specifically, the 
screw defect. However, the total chirality of the system however is not always the same for 
a given direction of the magnetic field, but changes from one sign of the chirality through 
zero to the other sign as the Fermi momentum or the magnitude of the magnetic flux is 
varied. In section 5, we write the (l+l)-dimensional effective field theory for these chiral 
states and couple it to an external electric field along the defect. This theory is anomalous 
as discussed above. This anomaly allows one to find the effective field theory in the bulk 
of the medium where the defect is embedded. It is given by Chern-Simons terms defined 
on the half-planes which have the line defect as the boundary. Finally, in the concluding 
section 6, we discuss other physical effects that can be deduced using this reformulation of 
defects in terms of non-trivial metric solutions of the gravitational equations. 



II. THE HAMILTONIAN ON THE PLANE 

In this section we will be discussing the motion of an electron in an external magnetic 
field along the z-direction. We will consider the Schrodinger equation for a spin half particle 
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with a magnetic moment interaction in this ambient space. 

The hne element describing the Euchdean geometry of is given by 

ds^ = {dr^ + r^de^) + dz'', (2.1) 

r,9, z being the cyhndrical coordinates. We will see later that various defect geometries can 
also be described by line elements like the above one. 
The spatial Laplacian is 

1 d d 1 (9^ 

^= [-^('"^) + ^^ + ^] • (2-2) 

We now split A according to 

A = At + Az = At + -^. (2.3) 

In the presence of a magnetic field along the 2:-direction and also a magnetic moment 
interaction, the electronic wavefunction on the plane satisfies the Pauli equation 

' ^ At - -a ■ B)^ = (2.4) 



' 2m* 2 

where At = DiDi is the transverse Laplacian constructed from the gauge- covariant deriva- 
tive 

A = di- iA, , (2.5) 

and and m* are the electronic magnetic moment and effective mass respectively. Notice 
that we have chosen our units by requiring the electronic charge e to be 1. We have also 
chosen a gauge where = 0. 

When = the Hamiltonian in (|2.4|) can be obtained by squaring the two-dimensional 



massless Dirac Hamiltonian. We will be dealing with this Dirac Hamiltonian below. 
For regular Cartesian coordinates our choice of the 7 matrices are 

7* = 7° = CT3; 7^^ = -ia2; 7^ = (2-6) 
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where ai are the standard Pauh matrices. From here one can readily find the 7 matrices for 
the polar basis : 



7° = 0^3, 7^ = (7 • r) 



( -e-^^^ 



\ 







7' = (7 • 0) 



J 



le-''^ 
ie'^ 



(2.7) 



We assume that B is time- independent. We can then choose a gauge where both Aq and 
are zero. In this gauge, the transverse Dirac Hamiltonian is identified via the equation 



2^\1> = where 



i^D = -^7"(7"^x + 7'^,) 

= 7°[-^(7 • 'r)dr - • 0){de - irAe)] 
f 



%e 



id 







le 



-ie 



l + tHm-rAe 







(2.8) 



This Hamiltonian is formally symmetric with respect to the measure r dr dO. 
Squaring this operator leads to the Hamiltonian ( |2.4| ) with ^J- = '■ 

1 



H 



2m' 



■HdHd. 



(2.9) 



Next we proceed to study the boundary conditions at the defect at r = appropriate 
for this Dirac operator. 

III. THE DIRAC EQUATION ON _ {g} 



The massless Dirac equation on a plane with a point removed admits a one-parameter 
family of boundary conditions at r = 0. This can be seen as follows. 

The Dirac operator —icr • defined (in the absence of the magnetic field ) on the 



two-dimensional plane, is self-adjoint on a domain T) if and only if |11 



rdr / d6 
Jo 



for ^ G P X G P, 



(3.1) 
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where 



Vt = r h d . 

dr r 89 



(3.2) 



Using Stokes' theorem and assuming that the wavef unctions fall off appropriately at spatial 
infinity, one ends up with the following equation at the boundary ( at r = 0): 



hm / de rx^A * = for * e X> x e P. 



(3.3) 



Here, 



A = —icr ' r 



( -ie''<^ ^ 



-ie'^ 



(3.4) 



The eigenvalues of the operator A are +i and —i and the corresponding eigenvectors nor- 
malized to 1 are 

/ . \ 

(3.5) 



1 

71 



^ 1 ^ 



and — ^ 
V2 



1 



respectively. It follows that the boundary conditions at r = compatible with self- 
adjointness of —icr • Vr are parametrized by e'^ e U{1) [ K being real] and are given 

by 



lim x/r*(r,6') = 



^1 



(0,9) ^a{9) 









( 1 


] 



























(3.6) 



where a{9) is any smooth function on S^. The parameter K (mod 27r) parametrizes the 
various boundary conditions. One can readily see that in terms of K, 



cos{^)MQ,9) - -te^^sm{^)MQ,9). 



(3.7) 



When X is or tt, this shows that the components ■02(0,^) or i/ji{0,9) vanish respectively 
and the solutions turn out to be "chiral" at the origin. 
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It is important that the Laplacian remains self-adjoint under the chosen conditions. 
We must therefore augment (|3.6| ) by the same types of boundary condition on the radial 
derivatives too, 



r->0 



^1^ 



pie 



+ e 



iK 



I 



\ 



1 ^ 



I 



(3.8) 



Here (Off) is any smooth function on S"^. [ It need not be a (6').] Equations ( |3.(j| ) and 
lead to 



/ 



rdr d9 



r-+0 



for ^ G P ^ X G P. (3.9) 



Here V = T>{e^^) is the set of functions fulfilling both the conditions (|3.6| ) and (|3.8| ). Thus 
the Laplacian is also self-adjoint with (|3.6| ) and (|378|) . 

We now include the magnetic field. It can easily be checked that the self-adjointness 
conditions discussed above remain unaffected as long as we are not dealing with singular 
magnetic field configurations. 

We will assume that our magnetic field is non-vanishing only over a compact region 
enclosing the defect. For simplicity we will in fact assume that the magnetic field is a 
constant within a circle of radius a and zero outside. ( The magnetic field can be represented 
by a delta function in the limit of a = 0, this case has been treated in detail by Moroz [Q.) 
Hence, our vector potential will be chosen to have the form 



A = Ardr + AgrdO 

Ar = 0, 



Br ~ $ ~ 

Ag = —9{a - r) + — 6'(r - a) 
I Zttt 



dA = B^rdr dO = B9{a - r)d^x = B9{a - r) rdrd9 



(3.10) 



where 



$ = Tia^B 



(3.11) 



is the total flux and 9 the step function. [ We denote it by 9 to avoid duplication of notation. 



As the magnetic field is cylindrically symmetric, the angular momentum for the Dirac 
particle is a good quantum number. Accordingly, the conditions ( p.6| ) and ( |3.8| ) can by 
analyzed in terms of the simultaneous eigenfunctions of Hd and the total angular momentum 
operator J = L + |cr3 where L = —i^ is the orbital angular momentum operator. The 
eigenfunction of the Dirac Hamiltonian Hd with total angular momentum j has the form 

(3.12) 



V 



where m = j — |. In terms of u^^^i and v^_^_l, the boundary condition ( |3.7[ ) reads as 

K K 



(3.13) 



The spectrum of angular momentum depends on the (quasi-)periodicity of the wave 
function under the 9 ^ 9 + 27r translation, 



^(r,9 + 27i) = e'^^(r. 



(3.14) 



which implies that 



A 

m G \- n, 

2tt 



n G 



The Hamiltonian Hd, restricted to the m-th partial wave \l/m(^ 
form 

/ 





H 



(m) 
D 



V 







D 




On squaring Hj^\ one obtains the Schrodinger operators 

DD^ ■ 



B, = -drirAe 
r 



(3.15) 
acquires the 



(3.16) 



(3.17) 



Aj^^ being the radial Laplacian including the gauge connection (which is a function of r 
only), 

(m — rAaY 



A. 



(m) 



(3.18) 



and Bz is the component of the magnetic field in the third direction, that is, in the direction 
normal to the plane. 

We now proceed to find zero energy solutions of the Dirac equation satisfying the bound- 
ary condition (|3.6|) . They automatically satisfy the boundary condition (|3.^ ) and hence are 
zero energy solutions of the Laplacian too. 

For a generic Ag^ the zero energy eigenfunctions of i/}^-* [ignoring the boundary condition 
( p.l3| ) and square integrability for the moment] are given by 



(3.19) 



There is a bound on the value of m from the requirement of square integrability of the 
eigenfunctions. Let us assume that $ is positive (i.e. the magnetic field is pointing along 
positive 2;-axis). Using the asymptotic forms 



Ae{r')dr' 







as r 0, 



(3.20) 



In r as r ^ oo 



which follow from ( p.lO| ), the solutions ( |3.19| ), the condition ( p.l3| ) and the requirements of 
square- integrability of the wavefunctions at r = and r — 5> oo, we obtain the wave function 





(3.21) 



and the following bounds on m : 



a) When e'^ ^ 1 



b) When e 



iK 



V " / 



1 $ 

— < m < ( — r 

2 

1 , $ 

— < m < ( 1 

2 - ^27r 



(3.22) 
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These bounds imply that 



$ > TT. 



(3.23) 



Thus, there are about — |] independent spin polarized zero energy solutions, [^] denoting 
the largest integer not exceeding ^. [ The precise number depends on the value of e*^ and 
the allowed values of m.] It is important to note that the parameter K determines whether 
the mode m = — | is a bound state or not. Of course this value of m is allowed only if that 
A = 7r(mod27r) as well. 

On the other hand, when $ is negative, we find the wave function 

/ . \ 







where m is now bounded as follows: 



c) When e'^ ^ -1 

d) When e^^ = -1 



1 $ 

-9 > > ^TT' 

Z ZTT 

1 $ 

— > m > — . 

2 - 27r 



As a consequence. 



(3.24) 



(3.25) 



$ < -TT. 



(3.26) 



Note that no zero energy bound states exist when $ < |7r|. 

We next proceed to extend this analysis to three dimensions and in particular to dislo- 
cations and screw defects. 

IV. DISCLINATION AND SCREW DISLOCATION 



A. Disclination 



The above analysis can be extended very easily to the case when we have a conical defect 
or a screw dislocation in our background. In the presence of a conical defect ( which one 
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can create by simply cutting out a wedge from a plane and identifying [pasting together] the 
edges or by inserting a wedge), the geometry of the defect can be represented by the metric 



ds' = {dr' + a'r'dO') + dz' (4.1) 

where the coordinates r, 6, z are the standard cylindrical coordinates [ with ^ = and 6 = 27r 
being identified as usual ] and 2n{l — a) is the opening angle of the cutout wedge. Note 
that < a < 1. In this coordinate system, the Laplacian can be written as 

1 <9 <9 1 92 
A = - — (r— ) H = — + — — . (4.2) 

r dr dr o?r'^ dff^ dz^ 

It is convenient to rescale the coordinates as follows : 

r — > i? = r, 
6) = a^, 

z^ Z = z. (4.3) 

The Laplacian then becomes 

1 (9 d 1 92 

Although A looks like the fiat space Laplacian, the range of is different, < < lua. 
As a consequence, the (quasi)periodicity of the wavefunction now reads 

^(i?,0 + 27rQ;) = e^^^(i?,0). (4.5) 

We also demand the identical periodicity condition on the radial derivative to ensure that 
A is self-adjoint. 

For the condition ( [4.5| ), the spectrum of the "angular momentum" gets quantized 

in units of -[m + ^\ , where m E Z^. The method of counting the bound states is the 

Q. V ' 271- / ' o 

same as in section 3, but for the change m — > ^. This incidentally also shows that the 
number of bound states can change in the presence of a disclination. 
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B. Screw Dislocation 



Next we move on to the screw dislocation. The screw dislocation can be characterized 
by the metric 



ds^ = {dr^ + a^^de^) + {dz + (3d9Y = gijdCd^^, 
d^^ = dr, di^ = de, di^ = dz 

where as usual r,6, z are the cylindrical coordinates for R^. 
The matrix [gij] and its inverse [g^^] can be written as 



[9^ 



1 

(«V + /32) (3 
(3 1 


















1 






^^2 2 




-/3 


(1+ f 



respectively. 



The Laplacian after a rescaling as in (4.3) is 
Ai = l^d,{g'^^d,) 



fir ArAi . Ar A _ ^ . 



^2 



'R^dR dR' R^^dQ adZ' dZ'^' 
g = det[gij\. 



(4.6) 



(4.7) 



(4.8) 



Wave functions in the domain of this Laplacian fulfill ( |3.7] ). 

It may be emphasized here that in the presence of a screw defect, the geodesies of the 
metric [gij] give the free propagation of low energy electrons.]^ 



^By making the coordinate change Z Z' = Z + ^Q, this problem can be mapped into that 
of a Laplacian on flat space, but with a non-trivial boundary condition. Specifically, the quasi- 
periodicity of the wavefunction under G ^ B + 27ra translates into 

^(i?, e, Z') = e'^^(i?, e + 27ra, Z' + 27r/3). 
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Next we add a magnetic field in the Z-direction. This is done as before, by replacing 
di by di — iAi where = 0; Aq = ^[^(a — R)^ + ^^d(R — a)] and also introducing a 
magnetic moment interaction. Note that the potential has been suitably scaled so that the 
total flux is still The Hamiltonian then reads 

1 



1 1 d d 1 2 d\ 



2m* 



-a- ■ B 



(4.9) 



where Ve = {de - iRAe - ^dz). 

The zero energy eigenfunctions of H are found as before, by first splitting H into a 
transverse part and a part along the Z-direction. The transverse Hamiltonian 



1 



1 



2m*'RdR"^"dR' R^ 
is still expressible as a square of a Dirac operator, 

^DDt ^ 

b^b 



2m* 



-a- • B 



(4.10) 



Hi 



H 



D 



\ 

b^ 



where 



2m 





2m* ^ 



(4.11) 



d 



ej 



(4.12) 



'dR R 

Requiring that D and be the adjoint of each other ( so that Hj^ is self-adjoint) leads 



once again to the same boundary conditions (3.6). 

We will make the ansatz that the wavefunction is of the form 



ikZ 



(4.13) 



Now, the zero energy modes of Hf) ( and hence those of Ht) [ignoring the boundary 
condition at the origin and square integrability for the moment] are given (up to an overall 
normalization factor) by 
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^^_(m+l)g/^''Ae(r')rfr'gi(.m+^fc+l)e 



(4.14) 



The corresponding functions for the full Hamiltonian ( [4 .91 ) are plane waves along the 
Z-direction and are given by 



£)/^-(m+l)g/jf Ae{r')*''gi(m+^A:+l)egifcZ 



/ 

(4.15) 



where m! = m + -k. 

a 

The case where both C and D are nonzero cannot occur as we shall see below. 

The quasi-periodicity condition (|4.5| ) leads to the following "quantization" conditions: 



i) When D = : am' = — + n, n E 

ii) When C = : a{m' + 1) = — + n, n e 

2tt 



(4.16) 



For positive $ , the requirement of square integrability at the origin (r 0) and at 
infinity (r oo) and satisfying the boundary condition ( |3.7| ) leads to the wavefunctions 





^(^,e,z) 



V 



(4.17) 



as well as to bounds analogous to (|3.22|) : 

a) When e''^ ^ 1 : < (m' - ^k) < - 1)), 

b) When e'^ = 1 : -- < (m' - ^k) < (— - 1). 

2 a 2n 

When $ is negative, the above requirements lead to the wavefunction 

' ^ 



(4.18) 



instead, and the bounds are now the analogs of ( [3.251 ): 



v 



(4.19) 
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c) When e'^ ^ -1 : -\ > im! - -k) > 

d) When e'^ = -1 : > (m' - -k) > -5-. (4.20) 

2 Q! 27r 

The energy eigenvalues associated with these wavefunctions are 

Eim,k)^- — . (4.21) 

The electronic states associated with the above wavefunctions will be loosely called zero 
mode electrons. The energy spectrum of these zero modes is given by 



k'^ a^{m' - mf (^ + 



n — am] 



2m* 2m*p^ 2m*/?2 ' ^ ' ^ 



where 

m for$>0(D = 0), 



m—\ 



(4.23) 

m + lfor$<0 (C = 0). 



There are degenerate levels which can be obtained by varying the integer n or alternatively 
by varying m' while holding n — am and hence m! — m fixed. 

In order to obtain the low energy effective action for the modes localized near the defect, 
we need to look at the behavior of the solutions near the defect. For $ > tt for example, 
the probability densities for the wave functions with m = — | are peaked around the origin 
and for those with m > —\ vanish at the origin. A similar situation prevails for $ < — tt. 
These peaked states are allowed only if e^^ = ±1. Apart from these particular states , 
the other states are localized at a finite distance away from the defect line. However, in a 
realistic situation , the defect line has a finite width and this will make some more states 
to be localized within the defect line, provided m is not too large and the magnetic field is 
large and can be treated as being uniform over the defect. This can be established from the 
Larmor formula, 

777 

r* = - (4.24) 

where r* is the radius of the Larmor orbit. The number of these low-energy excitations will 
be finite. 
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Though the above picture gives the single particle spectrum, one has to remember that 
in the many body picture for the fermions, one has to rather deal with the low energy 
excitations above Fermi surfaces. They can be approximated by linearizing the theory 
around each Fermi surface. 

Now we would like to see whether these excitations are chiral as a consequence of the 
chirality of the screw defect. In other words, we want to know if there can be some asymmetry 
between excitations "above" the two Fermi surfaces given by 



kp = ±\l2m*EF = ±\kF\, 



Ef = (4.25) 

The quasiparticle excitations are indeed chiral, as wc can sec in the following manner: 
We define their momenta as n = k — kp- For excitations "above" the Fermi surfaces ( which 
we will call particles ) one has, 

AE^ E-Ef^ —h > 0. (4.26) 

m* 

Thus for particle excitations at kp — \kF\-i one has k > ( upward or "right- movers" ) and 
for those "above" the Fermi surface at kp = —\kF\-, one has k < ( downward or "left 
movers" ). 

However, though the quasiparticles themselves are chiral, if the situation is such that the 
number of upward moving and downward moving quasiparticles are equal, there will not be 
any net chirality. Let us investigate this point further. 

For specificity, we assume that $ > tt and ^ > in the rest of the section. Similar 
considerations can be made for the remaining cases unless ^ = 0. We will comment on the 
latter case later. 

The relation between the linear momentum along the z-axis and the angular momentum 
is given by 

m^m' - -k (4.27) 
a 
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where m has to satisfy the bound ( |4.18| ). Also we have the fact that 



am' = — + n, ne (4.28) 
2tx 

This shows that for a fixed value of the angular momenta of two neighboring states differ 
by 

1 

Am' = -. (4.29) 

a 

In Figure |l| we have plotted the relation between k and m for different m's. 

The lowest energy electrons occupy the levels lying within the region bounded by the 
lines m = —\ and m = — 1 and k = ilkpl, which we will call the Fermi sea hereafter. 

We would like to determine the number of straight lines ( which are labeled by different 
values of m' ) crossing the two Fermi surfaces. These lines correspond to the low energy 
excitations of the systems above the Fermi surfaces. 

Note that the relation (|4.27| ) shows that as we increase k ( by applying external electric 



field, for instance ) , m decreases for fixed m'. As we have argued earlier that m is a measure 
of how far away from the defect the electron is localized , this shows that the electron is 
drawn into the defect as k is increased. ( Note that the direction of this current is dictated 
by the sign of £ and that this transport phenomenon will not occur if £ = ). 

Using relations ( [4. 271 ) and ( |4.29| ), one can now easily estimate the number of occupied 



states with a given momentum [the number of intersections of the graphs in Figure || with 
a horizontal line with a given k] , provided we know the maximum and minimum value of m 
for the states lying within our allowed region when k = 0, which we denote by rrimax and 
rrimin respectively. Thus, the number of excitation branches when A; = is 

A^o = a{mmax - rrimin) + 1- (4.30) 

Let us now define the following quantities 

1 

b = rrimax + 1 - {- 1) = {rrimax + 2 - — ). (4.31) 



The quantity a measures the horizontal distance between the hne m = — | and the 
occupied state which is just to the right of it for = 0, while h measures the horizontal 
distance between the line m = — 1) and the state which is just to the right of it for k = 0. 
(It moves into the Fermi sea when k is sufficiently increased). The value of a depends on 
the parameters A and a ands hence is fixed for a given material ( assuming that A cannot 
be varied ). On the other hand, b depends on the fiux $ and hence can be varied. Neither 
a nor b are functions oi kp- 

Now, as we increase k from zero to a positive value, the m values for the occupied states 
will change and those with m < — | will "fiow" out of the spectrum. A simple counting 
argument shows that the number of states that will "fiow" out will be given by the integer 



v. _ aa 
^ 

/3 



[Pk-aa]. (4.32) 



This equation can be obtained by noticing that the vertical spacing between two successive 
branches is Ak = ^ and that the minimum value of k required to move the left-most allowed 
state to the left of the bound m = — | is 

At the other extreme of the allowed region, there will states moving "into" the allowed 
bound. The number of these incoming states can be found by arguments similar to those 
above. This number is given by 

N2 = [pk-ab]. (4.33) 

Thus the number of branches appearing through the top Fermi surface at k = {kpl is 
given by 

Ntop = Nq- Ni + N2 = a{m.raax - rrimin) + 1 - [Plkpl - aa] + [(3\kF\ - ab]. (4.34) 

We can also evaluate the number of branches appearing through the bottom Fermi surface 
at k = —\kp\, employing procedure similar to above. However, this time we have to decrease 
the k value from zero to a negative value k < 0. Decreasing the k value will make the m value 
for various states to increase due to the relation ([4.27|) since ^ is positive by assumption. 
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Consequently, states will move away from the defect and some states with m < — | will 
move into the Fermi sea while some other states close to the bound m = — I will " flow" 
out of the Fermi sea. The number of branches appearing out of the Fermi surface k = —{kpl 
is then given by 

Nbottom = a(mmax " rUmin) + 1 + [/^I^fI + ao] - [/3\kF\ + ab]. (4.35) 

Thus the difference between the number of states from the top and bottom Fermi surfaces 
is 

77 = Ntop - Nbottom = [Plkpl - ab] + [Plkp] + ab] - ([/3|A;f| - aa] + [Plkpl + aa]). (4.36) 

This number t] determines the total chirality of the excitations. Note that when a = b , 
Nq = Ntop = Nbottom and 77 = 0. Thus, in this case the number of species remains the same 
irrespective of the value of {kpl- However, if a ^ b , r] can take any one of the three values 
, +1, and -1. This is shown in Figure ^ where 1] is plotted as a function of \kp\ for fixed 
values of aa and ab ( which we have chosen arbitrarily), r] can also change as the value of 
ab is changed ( by varying $ ) while {kpl and aa are kept fixed. This is shown in Figure |^. 
So, the direction of the current due to the quasiparticles will depend on the values of {kp], 
ab and aa. It depends in particular on $. 

We finally comment on the possibility f = 0. Since, |a| < 1, we then have /3 = which 
implies the absence of the screw defect. In this case, the straight lines in Figure |I| are vertical 
and the net chirality rj is always zero. 

The zero modes dominate the contribution to the effective action for the modes localized 
near the defect. We will analyze the low energy effective action describing these chiral mode 
further in Section V. 

V. TRANSPORT ALONG SOLENOIDAL DEFECTS AND ANOMALIES 

We have considered a situation where there are localized states in the presence of a 
solenoid with flux $. We will now consider the response of these trapped states to a con- 
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stant electric field along the axis of the solenoid ( that is, the defect). The low lying 
excitations about the Fermi surface corresponding to these states are described by massless 
chiral fermions localized on the defect. We proceed to discuss how their coupling to electro- 
magnetism leads to anomalous conduction properties. To set up the relevant background for 
this discussion, we first recall the reasons for the U{1) anomaly of chiral (1 + 1) dimensional 
gauged fermionic systems. 

Consider (1 + 1) dimensional chiral , say right-handed fermions, with field, i/jr = |(1 + 
■j5)ip. The action of this fermion field coupled to a constant electric field E{> 0) is given by 

S = I dtdxi^n^i^d^ + A^)iljR. (5.1) 

We choose the following 7 matrix convention: 

7* = ai; 7^' = -ia2, (5.2) 

so that 75 = 7*7^ = (73. The Dirac Hamiltonian for the 1+1 dimensional massless fermion is 

H = as{p-A,)-Ao. (5.3) 

Heisenberg's equation of motion gives, for the velocity v, 

V = j^{p - A,) = i[H,p - Ai] + -^^{p - A,) = -E. (5.4) 

The derivation of the chiral U{1) anomaly for the action (|5.1| ) by Nielsen and Ninomiya 



14 1 is as follows. The density of states p{v) for a massless (1 + 1) dimensional particle is 



controlled by the following formula for its variation : 



Sp{v) = —Sv. (5.5) 

271" 

Hence the production rate of these excitations per unit length is 

Nn(t) = I = J-i = f (5.6) 

dt In In 

where Nnit) is the number of fermions per unit length and we have used ( |5.4| ). As these 
excitations carry charge, there is a non- conservation of electric charge Qnit) carried by the 
right-handed fermions in the (1+1) dimensional theory, the rate of change of QR{t) being 
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Qnit) = eNR{t)L = -^L (5.7) 

Zn 

where L is the size of the sample. In this way, we obtain the U{1) anomaly. 

One can show in a similar way that for a left-handed chiral fermion, the rate of change 
of its charge is given by 

QL{t) = eNL{t)L = (5.8) 

where Nilt) is the number density of the left-handed fermions. 

Thus if there are rii species of right-handed fermions and n2 species of left-handed 
fermions, then the rate of change of electric charge is given by 

^AQnii) + Qdt)) = Qtotai = -(ni - n2)^L. (5.9) 

at ZTT 

As we are going to treat the different branches of the quasiparticle excitations as different 
species for the screw defect situation, one can readily see that 

{rii - n2) = ry (5.10) 

where t] is given by ( [4.36| ). 

But one knows that the system as a whole has charge conservation. Therefore, when 
?7 7^ 0, it must be that charge from outside flows into or away from the defect, and there 
must be non-conservation of charge outside the defect too. This fact can be encoded in an 
anomalous Chern-Simons "effective" action, representing the electronic degrees of freedom 
that are not localized on the defect. The connection between the anomaly (|5.7|) and the 



Chern-Simons term is elegantly demonstrated in |0( See also |15|). 

Using their results, it can be shown that the contribution appearing on the right-hand 
side of the equation (|5.7| ) is cancelled by the bulk action given by 

V 



S, 



bulk 



dt rdr dzAdA 



(5.11) 



where rj is given by the equation ( |4.3(j| ). The Chern-Simons three- form is being integrated 



here on an infinite half-plane which has the defect as the boundary, as shown in Figure 
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The gauge variation of the Chern-Simons term produces a surface term which is exactly 
cancelled by the variation of the chiral anomaly on the defect PJT5[1, the coefficient in 
( p.llp is chosen that this cancellation occurs. 

The fact that the Chern-Simons coefficient is determined by the anomaly fixes of the 
transverse conductivity of the defect, there being an exchange of current with the bulk when 
an electric field is applied along the line defect. This transverse conductivity is given by 
X 271 = the 2n coming from the (p integration. [ The first factor comes from varying 



A [|^. Recall that we have set the electric charge e = 1, if we had not done so, this 

2 

conductivity would have been As the value of rj can change by ±1 only, the system 
will behave somewhat like the integer Quantum Hall samples thought here the conductivity 
is bounded and the sign of the current is fixed by rj . 



VI. CONCLUSION 

In this paper we have studied the effect of line defects on the transport properties of 
solids in the presence of a magnetic fiux parallel to the line. This was done by describing 
the defects in a crystal in terms of a non-trivial metric and then studying the propagation of 
electrons in this background metric. We have also discussed the low energy effective action 
for this system. It is a 1 + 1 dimensional chiral fermionic action localized on the defect. 
The coupling of this system to an external electric field parallel to the defect leads to an 
anomaly and induces a bulk action involving the Chern-Simons term. The resultant total 
action in turn leads a phenomenon similar to the transverse Hall effect with a quantized 
conductivity. However, the sign of the conductivity, which is proportional to rj, undergoes 
flips as one changes the value of the magnetic flux or the Fermi momentum . 

There are numerous interesting ideas suggested by the description of defects in terms of 
metrics and we will now discuss a few of them. It has been shown that in the background 
of disclination defects, electric charges experience a non-trivial force which arises due to the 
conical defect angle and is directed towards, or away from, the defect line. For a positive 
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conical angle ( where a wedge has been taken out ), this force turns out to be attractive ( 
that is, towards the defect line ) and leads to bound states of electrons and defects [Q]. It 
would be interesting to see how this force is modified by the screw dislocation. 

We have assumed here that the magnetic field is localized in the neighborhood of the 
defect. In a type-II superconductor such a scenario can actually happen. The strain energy 
due to the defect will force the region around the defect to be in the normal state. Then it 
is favorable energetically for the magnetic flux tubes to be trapped within the screw defects 
( if they are present in the material) [l^. This is basically due to Meissner effect. Our 
results, therefore, would be interesting in the study of fiux-line trapping in superconductors 



by defects [jT6|. 



In this paper we have discussed the fact that electronic conductivity transverse to the 
defect is fixed in quantized units because the low energy effective theory of the electrons 
localized on the defects has an electromagnetic anomaly. It is known however that there is 
a gravitational anomaly in the same (1+1) dimensional system. It will lead to a transverse 
"gravitational Hall effect" , and this will manifest itself in anomalous elastic and vibrational 
properties. Work is in progress to analyze in detail how this affects the properties of the 
system. 

The fact that static defects can be described in terms of stationary solutions to 2+1 
dimensional Einstein action for gravity with matter gives rise to the interesting possibility 
that the dynamics of these defects can be modeled in terms of a fully dynamical theory 
of gravity. It would be interesting to carry out the study of dynamics of defects in the 
language of quantum gravity in two and three dimensions. In particular it suggests the 
interesting possibility that systems with defects could provide us with analog simulations of 
many situations in quantum gravity that are now being studied using elaborate computer 
simulations. 

The interplay of ideas from two widely different fields like condensed matter physics and 
quantum gravity can lead to lots of new results of the kind we have just discussed. It would 
be very interesting to use these condensed matter systems as low energy experimental probes 
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into what is conventionally regarded as the domain of Planck scale physics. Similar analog 
probes of cosmological defects have already led to rewarding results ( |p!7| ) 
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FIG. 2. Plot of r] vs. /3\kF\ for aa = .23 and ab 
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FIG. 3. Plot of r] vs. $ with aa = .37 and l3\kF\ = 5.7 and m. 
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The defect line 




FIG. 4. Half-planes with the defect as the boundary. 
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